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WHISTLER WAVES OBSERVED UPSTREAM FROM 
COLLI SONLESS SHOCKS 


D.H. Fairfield 


May 1973 



Abstract 


Waves in the frequency range 0.5 - 4 hz have been studied 
in the region upstream of the earths bow shock using data from 
the fluxgate magnetic field experiment on IMP-6. Such waves 
are invariably detected adjacent to the shock, persisting 
upstream for intervals often less than a minute but occasionally 
of the order of many hours. Analysis of 150 examples of these 
waves during a three month interval indicates that amplitudes 
are generally less than 1 or 2 gammas (AB/B - 0.2) and propaga- 
tion directions generally make angles of between 20° and 40° 
with the field direction. The waves as measured in the space- 
craft frame of reference are either left or right hand 
polarized with respect to the average field direction. The 
left handed waves generally have lower frequencies than the 
right handed waves and the left handed frequencies never 
exceed 2. 5 hz. The measured sense of polarization is found to 
depend on the propagation direction (or alternatively the 
field direction) relative to the solar wind direction. 
Propagation-solar wind angles greater than approximately 40° 

(or field angles greater than 55°) correspond to right hand 
polarization and smaller angles correspond to left handed 
polarization. It is concluded that the observed waves are 



right handed waves in the plasma frame of reference with wave- 
lengths of approximately 100 km propagating upstream in the 
whistler mode. Doppler shifting reduces the observed 
frequencies in the spacecraft frame and reverses the ob- 
served polarization for those waves propagating more directly 
upstream. Similar waves are seen ahead of most interplanetary 
shocks . 
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I. Introduction 


Large amplitude quasi-sinusoidal waves upstream from the 
earth's bow shock appear to fall naturally into two frequency 
ranges of approximatelv._. 01-.05 Hz and 0. 5-4. hz_.__ The lower 
frequency waves have been studied by Fairfield (1969) and 
Greenstadt et al . (1970b) and their properties are fairly 
well-known. They are primarily transverse waves of several 
gammas amplitude and they exhibit left hand polarization 
in the spacecraft frame of reference. The waves are associated 
with the bow shock since they are only observed on inter- 
planetary field lines that connect with the shock (Fairfield, 
1969). Furthermore these waves appear to be instrumental in 
determining the observed character of the bow shock (Greenstadt 
et al . , 1970a; 1970c; Greenstadt, 1972). It has been suggested 
that these waves are locally generated in the solar 
wind by protons moving upstream from the shock (Fairfield, 

1969; Barnes 1970; Greenstadt et al . 1970b). 

Characteristics of the 0.5 — 4. hz waves are not as clearly 
known. Heppner et al . (1967) first noted the presence of such 
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waves which were found to occur within a very few minutes 
of a bow shock crossing. These authors distinguished between 
frequency regions of .5 to 2.5 hz and those greater than 3 hz 
but were unable to determine the polarization of these up- 
stream waves. Propagation from the bow shock in the whistler 
mode was proposed as the likely explanation for these waves, 
primarily because this is the only mode with high enough group 
velocity to propagate upstream against the solar wind which 
streams into the bow shock at approximately 400 km/sec. 

Holzer et al . (1972) studied three bow shock crossings 

which occurred within a time interval of eight minutes and 
found that in all three cases the high frequency type waves 
exhibited right hand polarization. This observation demonstrated 
that the waves were not standing with respect to the bow shock 
(Tidman and Northrop , 1968) since a reversal of the direction of 
the shock motion did not produce a reversal in the observed 
polarization . 

Russell et al . (1971) studied wave packets with frequencies 

near 0.4 hz which occur upstream of and are associated with the 
bow shock but which are not restricted to the immediate 
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vicinity of the shock. These authors also noted the 
occasional presence of longer wave trains with frequencies 
near those of the wave packets but they did not include these 
waves in their study. The packets were generally observed to 
exhibit left hand polarization, but it was argued that 
the polarization was right handed in the solar wind frame 
of reference with the observed reversal in polarization 
being due to movement of this reference frame past the space- 
craft at the solar wind velocity. The authors rejected the 
explanation of whistler propagation on the basis of slightly 
noncircular (i.e. elliptical) polarization of the waves. 
Alternatively they proposed a particle origin similar to 
that proposed for the low frequency waves. Subsequently Wu 
(1972) has proposed an echo phenomenon to explain the ob- 
servations of Russell et al . and Hasegawa (1972) has proposed 
a whistler instability. 

An example illustrating the presence of both .01 - .05 hz 
and 0.5 - 4. hz . frequency waves is presented in Figure 1. 

In the remainder of this paper the terms ’’high" and "low" 
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frequency waves will be used to refer to these two classes 
of waves. In the top panel 30 minutes of data is shown 
which was taken during: an interval ending- 90 minutes prior 
to the subs equent_pbservat ion of the bow shock at a position 
more than 1.5 R downstream. The field magnitude F and solar 

JCi 

ecliptic latitude and longitude angles 9 and 0 are plotted 
for field averages spaced 15.36 seconds apart. Also, the 
square root of the sum of the squares of the three component 
standard deviations associated with the average are plotted 
as the quantity 6 . In the bottom panel individual measure- 
ments for the two minute interval from 19h 44m to 19h 46m are 
shown in the same format. The top panel clearly shows the 
presence of low frequency waves with periods near 50 seconds 
and the bottom panel illustrates the simultaneous presence 
of ~1 second period waves. Although both high and low frequency 
waves are simultanequly present in figure 1, either type may be 
seen in the absence of the other. This fact along with the 
existance of two characteristic freauencv ranges suggests the 
possibility that high and low frequency waves are two distinct 
phenomena with separate origins. 
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The present paper reports on the further investiga- 
tion of the higher frequency waves, particularly as they 
occur adjacent to the bow shock. It will be argued that 
a systematic analysis of the observations of both left and 
right hand polarizations confirms that they are indeed 
whistler mode waves. 

I I . Experime nt 

The IMP 6 (Explorer 43) spacecraft was launched on 
March 13, 1971 into an eccentric orbit with the initial apogee 
at a geocentric distance of 33.12 R (earth radii) at a 

£1 

longitude 15° east of the solar direction. The inclination of the 
orbit was 28.7° and the period was 4.18 days. During the 
first 30 orbits bow shock crossings occurred in the northern 
hemisphere within 12 R of the ecliptic plane between the 
solar ecliptic longitudes of 250° and 35°. The spacecraft 
was spin stabilized with spin axis approximately perpendicular 
to the ecliptic plane. The spin period was 11.1 seconds 
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during the first six months of the lifetime. The spacecraft 


was designed to operate in 400 and 1600 bit per second tele- 
metry modes but since it has operated virtually exclusively 
in the higher mode during the first two years of its lifetime, 
sampling frequencies cited below will refer to the higher 
mode. 

The IMP 6 magnetometer experiment (Seek et al , 1973) con- 
sisted of a triaxial fluxgate magnetometer with one sensor 
parallel and two sensors perpendicular to the spacecraft 
spin axis. The sensor were mounted at the end of a 13 foot 
boom to minimize the effects of spacecraft magnetic fields. 
Prelaunch testing indicated that any spacecraft field was less 
than 0.1 Y at the sensor position (Harris, 1972). 

The experiment operated in any one of the four ranges 
+16, +48, +144 and +432y. The three orthogonal sensors were 
sampled within an interval of 2.0 milliseconds and successive 
vector samples were taken at 80 millisecond intervals. The 
analog outputs from each sensor were digitized onboard the 
spacecraft with an 8 bit A-D converter. Full 8 bit words 
were telemetered to the ground every 16 measurements (1.28 sec) 
but for the intervening measurements differences were taken 
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between successive samples and only these differences were 
supplied to the telemetry stream. These differences were 
restricted to 4 bit words, thus allowing a 78% increase in 
the sampling rate for the same telemetry assignment. Summing 
these differences during ground processing allowed reconstruc- 
tion of the full bit words. In addition, the use of 16 
differences introduced redundancy since the sum of 16 
differences should be identical to the following full bit 
word. Rare instances when this sum did not yield the next 
full bit word indicated an inconsistency due either to telemetry 
noise or to exceeding the capacity of the 4 bit difference. 
Instances of exceeding the capacity of the four bit word 
were minimized by automatically switching the instrument to 
a less sensitive range when this capacity was exceeded. 

The instrument also automatically switched ranges according 
to criteria based on field magnitude. Return to a more 
sensitive range required both low enough fields and small 
enough differences (Seek et al . , 1973). Quantization step size 
due to digitization in the four ranges was + . 06,+ . 19 ,+ . 56 +.69 y< 



Zero levels of the sensors were determined by utilizing 


the spin of the spacecraft to reverse the orientation of the 

sensors perpendicular to the spin aixs. Every 46.5 hours 

o 

the sensor set automatically "flipped" 90 about the axis de- 
fined by a perpendicular sensor oriented along the radial 
direction. This action reversed the relative positions of 
two sensors parallel and perpendicular to the spin axis and 
allowed a determination of the zero level of the third sensor 
during the subsequent 46.5 hours. These zero level determina- 
tions were carried out independently for each range. The 
frequency of the zero determinations depended on the quietness 
of the ambient field and the expected time the sensor would 
be in a given range. Individual zero determinations were made 
at intervals which varied from every minute to every 10 minutes, 
but all data from half an orbit were combined and updated in 
the data processing program on a twice per orbit basis. The 
zero levels varied by less than 2y from their launch values 
during the first 16 months of operation. The useful accuracies 
of the zero levels after they are corrected by the above pro- 
cedures are estimated as + O.ly in the lowest range and + 0.3Y 
in the highest range. 

Every 11.7 hours a calibration field of 5y (50y in the 
higher two ranges) was added along each sensor axis to check 
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the sensitivity of the magnetometer. No change in the pre- 
launch values have been necessary during the first 1.5 years 
of operation. 

The bandpass of the instrument (3 db attenuation at 
7 Hz) was choosen to correspond to the 6.25 hz Nyquist 
frequency of the experiment. This insures that the 0.5 - 4. Hz 
waves reported in this paper are being observed at their 
true frequency and that aliasing is not important. 

Ill . Analysis and Results 

Observations . The present study considered approximately 
150 bow shock crossings which occurred on 19 of the first 30 
orbits of IMP 6. These cases were typical high mach number 
shocks with rare low mach number shocks (Fairfield, 1971; Formisano 
et al,1971) wifh(M^ ~ 3) being eliminated from consideration. 

In most cases the bow shock is readily identified by an increase 
in the magnetic field strength by a factor of approximately 3. 

In cases when the interplanetary field tends to be aligned 
with the shock normal, low frequency upstream waves are present 
(Fairfield, 1969; Greenstadt, 1970b) and the characteristic 
shock increase often is obscured by large amplitude waves or 
pulses (Greenstadt et al., 1970a, 1970b; Greenstadt, 1972). 

Such data sometimes occur continually or recurrently 
over an interval of several tens of minutes and 
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the determination of the precise moment of "a shock crossing" 
becomes rather subjective. The present study has utilized data 
sampled at .08 sec intervals and the policy has been to consider 
as individual shock crossings those increases that are clearly 
separated by low magnitude interplanetary data for at least 
30 seconds. 

The set of shock crossings was divided into two groups 
of "clean" and "pulation" crossings according to the presence 
or absence of low frequency upstream waves and (usually 
equivalently) the presence or absence of a ragged and pulsating 
character of the data. Examples of the two types of crossings 
are shown in Figure 2. The example from April 20 shows a clean 
shock crossing occurring at Ilh27m28s where no low frequency upstream 
waves are present. The example from April 12 shows a pulsation 
crossing which would probably be located at 5h40m25s, however 
the additional low field region at 5h40ml2s and other variations 
at earlier times (not shown) makes the identification somewhat 
ambiguous. The presence of irregular low frequency upstream 
waves (periods approximately 20 sec) can also be seen in 
this example. 
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For all of the shocks considered, the angle between 
the field (generally averaged over a 30 or 60 second interval) 
and the shock normal was calculated. The shock normal was 
computed from the average two dimensional bow shock determined 
by Fairfield (1971) in the three dimensional form given by 
Scudder et al . (1973). Figure 3 illustrates the histograms 
of the pulsation and clean crossing as a function of this angle 
between the field and the normal. The clean crossings are 
restricted almost exclusively to angles greater than 50°, 
whereas the pulsation crossings are primarily confined to 
angles less than this value. This result is consistent with 
the results of Fairfield (1969) and corresponds closely to 
the results of Greenstadt (1970b) when presented in terms of 
an average angle (Northrop and Birmingham, 1973) . The lack 
of pulsation cases for low angles is due to the decreased, 
probability of the field falling within a small solid angle 
surrounding the normal direction, as well as the fact that a 
field aligned with the shock normal would not be expected to 
increase at the shock and produce the characteristic signature 
necessary to identify the shock with only magnetic field data. 
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Figure 2 also serves to illustrate the higher frequency 
waves which are the primary subject of this paper. In each 
of the examples of Figure 2 these waves are present throughout 
the upstream region though the frequency is clearly higher on 
April 20 than on April 12. 

In many cases, however, the wave amplitude decreases 
as the time from the shock crossing increases. This 
time duration, outside of which the wave amplitude is imperceptable 
on the plots ,is typically less than a few minutes for the clean 
crossings, although occasionally the duration is of the order 
of hours or even tens of hours. For the pulsating crossings 
the frequent disappearance and reappearance of the waves 
makes it much more difficult to specify a time duration for 
the waves. To quantify this analysis the time duration adjacent 
to the shock crossing during which the high frequency waves 
persisted was scaled from each of the clean shock crossings. 

Since this time is presumably related to some damping length, 
cases of multiple crossings (within a few minutes) where 
the spacecraft remained in the vicinity of the shock were 
omitted from this analysis. In figure 4 this persistance 
time, At, is plotted for the clear shocks versus the angle 
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between the field and the shock normal. Although points do 
not appear at low angles due to the omission of the pulsating 
shocks, there is still a clear tendency for the waves to 
persist for a shorter duration upstream when the field- normal 
angle is large. This persistance time cannot be converted to 
a distance from the shock since the relative velocity 
between the spacecraft and the moving shock is not known. 

To determine the frequencies of the waves adjacent to 
and upstream of the shock, power spectra were calculated for 
130 shocks where the waves were present for an adequate amount 
of time. The spectra utilized one minute of data when 
available (30 seconds when necessary) and had 25 degrees of 
freedom. The analysis was carried out in a coordinate system 
where Z was the average field direction and X and Y were 
transverse to this direction. In addition to the four single 
component spectra (three components and the magnitude) the 
coherence and phase between two transverse components was 
calculated as a function of frequency (Bendat and Piersol , 
1966) . The coherence was invariably high near the peak in 
the spectrum and the corresponding phase was invariably near 
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90° or 270°. These phases indicate either left or right hand 
polarization respectively as the sense of rotation of the wave 
vector relative to the average field direction. 

Power spectra corresponding to the examples of Figure 2 

are shown in Figure 5 for a field component transverse to the 

average field direction. The dashed line spectrum (April 12) 

exhibits a rapid decrease for frequencies above the peak. 

Such left handed spectra often reach the digitization noise 
-4 2 

level of 2 x 10 y /hz for frequencies less than 6 hz. The 
solid line spectrum corresponding to the right hand waves is 
typical in that the peak is at higher frequencies than that 
for left hand waves. Usually the decrease in power for right 
hand waves is less abrupt and it seldom reaches the digitization 
noise level of the left handed cases. In most right hand 
cases there is a peak followed by a decreasing spectrum which 
suggests that the important part of the spectrum is covered by 
the instrumental range. This observation is supported by 
the work of Olson et al (1969) who found that the spectrum gener- 
ally falls rapidly above approximately 5 or 10 hz. On rare occasions 
when the polarization is right handed the present work reveals 
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only a minimal decrease at high frequencies (due partly to 
instrument attenuation) and in these rare cases there may be 
a significant spectral enhancement beyond 6.25 Hz. The spectra 
of left hand polarized waves always decreases sharply and 
there is no suggestion that left handed waves have ever fallen 
outside of the instrumental range. 

Of the 130 spectra calculated adjacent to the bow 
shock some showed a peak which was even more pronounced than those 
of Figure 5, but many others exhibited peaks which were broader 
and less well defined than Figure 5. In the great majority of cases, 
however, either left handed or right handed polarization could 
be identified as predominating. In a few cases there was 
evidence for left handed polarization at low frequencies and 
right handed polarization for higher frequencies. 

Altogether 112 spectra were processed where a polariza- 
tion and at least an approximate frequency could be determined 
for the upstream waves adjacent to the shock. Figure 6 
illustrates the peak frequency of the waves plotted versus the 
angle between the field and the X axis (the approximate direc- 
tion of the solar wind) . Right hand polarizations have been 
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plotted as positive frequencies and left hand polarizations 
as negative frequencies. The lack of points with ]f| < .3 hz 


is due to 

the 

fact 

that 

the spectra computed 

in 

this study 

did 

not 

have 

adequate resolution in 

this 


low frequency regime. Furthermore, frequencies below 0.3 hz 
begin to overlap the .01 - .05 hz frequency range where most 
of the observed waves apparently have a different origin. 

Clearly an angle of 55° + 5° in figure 6 separates the cases 
of right and left polarization. Each group of polarization 
in figure 6 contains many cases of interplanetary to magneto- 
sheath traversals as well as magnetosheath to interplanetary 
traversals. Since the interplanetary field orientation determines 
the sense of polarization so well and since the field orientation 
cannot be well correlated with the direction of shock motion, this 
figure can be cited as confirmation of the result of Holzer et al . 
(1972) that the polarization does not depend on the relative direc- 
tion of motion between the shock and the satellite. 

Another characteristic of prime importance in the study 
of waves is their propagation vector k. This direction can 
be determined by a technique originally developed by Sonnerup 
and Cahill (1967) for studying the magnetopause, but 
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subsequently applied to waves by many authors (e.g. Siscoe 
et al . , 1967; Russell et al . , 1971; Holzer et al . , 1972). 

The procedure is to diagonalize the variance matrix in order 

to determine the direction of minimum variance which is 

a , 

assumed to be the£k direction of the waves. The procedure 

gives no information about the magnitude of k and also leaves 

-3* 

an ambiguity in the sign of k which is surpressed by assuming 
the propagation is in the upstream direction. This variance 
analysis was carried out on the' upstream wave events for which 
spectra were available. In each case the technique was applied to 
successive intervals of 4 seconds (several cycles of the 
wave) covering an interval of typically half a minute. This 

A : : 

procedure allowed intercomparison of the computed k vectors from the 
various segments. In 84 cases the individually determined 

A 

k vectors were judged to agree well enough among themselves 
(i.e.to fall roughly within a cone of 15° half angle) to warrant 
using the event. In these acceptable cases both the angles 
between the field and the V vector and the k vector and the 
X axis were averaged. 

Figure 7 is a histogram showing the number of cases with 
different values of the angle between the upstream field vector 
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and &. There is a strong preference for waves to propagate 
at an angle between 20° and 40° to the field with a clear 
tendency for the waves to neither propagate along nor 
perpendicular to the field. The fact that the waves generally 
propagate preferentially at an angle within 40° of the field 

A 

direction suggests that the k vector might possibly be 
substituted for F in ordering the data of Figure 6. 

Figure 8 is the figure similar to Figure 6 only with 
the abscissa now being the angle between Is and the approximate 
solar wind direction, x. The points represent a subset of those 

A 

of Figure 6 corresponding to those cases where a k vector 
could be determined. The crosses represent occurrences 
of the waves at locations more than 2 hours upstream of the 
nearest shock crossing. Indeed the k-x angle orders the data 
about as well as the F-x angle of Figure 6, particularly 
considering the greater uncertainties in the determination 
of £. Now a kx angle of about 40° separates the regions of 
positive frequency (right hand polarization) and negative 
frequency (left hand polarization) . It also appears that 
right hand waves occur at frequencies above 2.5 hz, whereas 
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left-hand waves do not (do not fall below -2.5 hz in Figure 
8) . It will be argued in the next section that Figure 8 is the 
key to confirming the origin of the waves as whistler mode 
propagation from the bow shock. 

It should be pointed out that due to the geometry of the 
bow shock there is a tendency for the shock normals to point 
in a generally upstream direction. This tendency means that 
there is a strong statistical tendency toward the condition 
for low frequency upstream waves and pulsating shocks (the 
field aligned generally with n: see Figure 3) to occur at 

. , J 

the same time as the condition for left handed high frequency 
waves (the field generally aligned with x: see Figure 6) . 

This statistical correlation can be violated, however, and 
pulsation shocks can accompany right handed waves and clean 
shocks can accompany left handed waves. This latter case is 
illustrated by Figure 9 which shows the detailed data for a 
shock crossing on April 3, 1971. These data appear to be 
typical of clean shock crossings but they are unusual in that 
the upstream waves are left hand polarized rather than right 
hand polarized as they are for most clean shock crossings. 
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By looking at the various angles we can understand this 
polarization (see insert in Figure 9) . The field makes 
an angle of 50° with the shock normal which is apparently 
just large enough to exclude the low frequency waves and 
create the "clean" shock. At the same time the field 
(which happens to be oriented almost perpendicular to the 
average spiral angle) makes an angle of only 37° with the 
X axis (kx angle is 12°) and therefore the waves are left 
handed. It appears that exceptions to the general correla- 
tion of pulsating crossings with left handed high frequency 
waves and clean crossings with right handed high frequency 
waves can be explained in this manner. There appears to 
be no reason to suggest that the statistical correlation 
has any independent significance. 

To illustrate that high frequency waves are related to 
shocks other than the earth's bow shock, Figure 10 presents the 

data at the time of two interplanetary shocks. On May 30 
a very clear wave train is apparent ahead of the shock and 
power spectra indicate a clear peak at 1.9 Hz with right 
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handed polarization. On May 17, higher frequency (4.2 Hz) 

right handed waves are present upstream from the shock. In 

a o o 

these two cases k makes angles of 46 and 65 respectively 

with the solar wind direction. At the time of three of four other 

available interplanetary shock waves, similar waves are 

seen, although in two cases the amplitude is too small and 

the duration (~4 seconds) too short to permit further analysis. 

The remaining analyzable case had a frequency of 1.0 Hz, 

exhibited right hand polarization, and had a kx angle of 

73°. Only on one interplanetary shock with the field aligned 

near the shock normal were no waves observed. 

Interpretation . The association of 0.5 - 4. hz waves 
with the bow shock and their apparent damping with distance 
from it suggests that the waves may be propagating away from 
the shock. As was pointed out by Heppner et al . (1967), the 
whistler mode is an obvious candidate for their explanation 
since it has a group velocity greater than the solar wind 
velocity. The wave frequencies under consideration are well 
below the electron and ion plasma frequencies and the electron 
gyrofrequency but above the ion gyrof requency so the simple 
cold plasma dispersion relation(e. g . Stix 1962) for right 
hand polarized whistlers becomes 
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where m ' is the wave frequency in the plasma frame, ou is 
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the electron plasma frequency (co = — ) , o is the 
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electron gyrofrequency (Q = ) , 0 is the propagation angle 
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with respect to the field direction (the angle in Figure 7), 
n is the plasma density and k = 2tt/\ . 

Waves observed at the spacecraft will have undergone a 
doppler shift due to motion of their frame of reference past 
the spacecraft. They will be detected at frequency tu where 


u> = w ' + k* (2) 

and V is the solar wind velocity. With the assumption 
that propagation is in the upstream direction, we may re- 
write 2 as 


sw 

f = f' - cos (kx) (2*) 

\ 

where the kx angle is that of Figure 8. 

The observed frequency f will be the plasma frame f’ reduced by 
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a factor depending on the solar wind velocity, the wavelength 
and the direction of propagation. Physically a negative f 
in equation 2 ' means that the wave whose perturbation vector 
traces a right handed helix as it propagates upstream in the 
plasma frame is being convected downstream past the spacecraft 
in such a way that the observed polarization is reversed. For 
this reversal to take place the phase velocity of the wave 
must be less than the solar wind velocity but at the same 
time the group velocity must be greater than the solar wind 
velocity to enable the waves to propagate away from the shock. 
(Actually the above is strictly true only at the subsolar 
point for propagation in an exactly upstream direction. The 
problem should actually be approached three dimensionally 
in which case the shock normal is introduced into the 
problem and the component of the group velocity along the 
normal is required to exceed the solar wind velocity component 
along the normal. Only the one dimensional argument will be 
presented here, but it should suffice for the order of 
magnitude arguments given below. ) 
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The solid curve in Figure 11 is a plot of the dispersion 


relation of equation 1 (f in hertz vs k) for a typical set of 
solar wind parameters. In addition the phase velocity 

(JJ d (i) 

(V _ —'l and group velocities (V = — ) have been calculated 
ph _ k g dk 

from the dispersion relation and are shown as dashed lines. 


These velocities correspond to the scale on the right. It 

is of interest to examine the phase and group velocities relative 

to a typical solar wind velocity of 400 km/sec to see whether 

the conditions of V > 400 km/sec, V , <- 400 km/sec can be met 

S ph 

simultaneously. For k>.65 (\<10 km;f , ’>40 hz) the group velocity is 
seen to be less than 400 km/sec and waves could not even propagate 


away from the shock. For . 09 < k < .65 (10 < \ < 70 km; 

5.7 < f 7 < 40 hz) the waves could propagate upstream, but 
their phase velocity is also larger than the solar wind velocity 
and consequently no left hand polarizations would be observed. 
For k <- .043 (\ > 148 km; f ' < 1.4 hz) the group velocity 
is again too low to allow upstream propagation. (Care must 
be exercized at very low frequencies where the neglected ion 
term in the dispersion relation becomes important) . 
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Only in the very limited frequency range where .043 < k . 09 
(70 < > < 148 km; 1.4c f' < 5.7 hz'i would we expect to see 
both left and right handed polarization because only in this 
frequency range is the group velocity greater than 400 km/sec 
and the phase velocity less than 400 km/sec. This 
theoretical frequency range corresponds very closely to the 
frequency range where both left and right handed waves are 
observed. 



This prediction is exactly the result illustrated in Figure 
8 which shows that right handed waves are observed for large 
angles (small values of the cosine) and left handed waves for 

small angles (large cosines). 

On the basis of Figure 8 it can be said that zero 
frequency waves should occur at a kx angle of ~40° when the 


phase velocity of the wave equals the component of the solar 
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A 

wind velocity in the k direction. Assuming a solar wind 
velocity of 400 km/sec, we obtain an average phase velocity 
of 

(W o 

— = 400 cos 40 = 306 km/sec. 

k 

Referring to Figure 11, we find that this velocity corresponds 
to a frequency in the plasma frame of 3.2 hz which in turn, 
through the dispersion relation, corresponds to a wavelength 
of 95 km. 

It is interesting to note that this frequency of 3.2 hz 
is the lower hybrid frequency in an average field of 5y, 
raising the question of whether the lower hybrid frequency 
is a significant frequency. In principle, it is possible 
to see whether the waves always occur at the lower hybrid 
frequency by solving equations 1 and 2 simultaneously with 
the measured parameters from individual cases. In practice, 
this solution involves a cubic equation whose solution has 
not yet been attempted. Individual cases have been numerically 
investigated, however, (using the measured solar wind 
parameter obtained from by the Los Almos plasma experiment on 
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board the same spacecraft) and it appears that although 
individual determinations of the plasma frame frequencies 
would cluster about the lower hybrid frequency, they do not 
correspond to it on a case by case basis. 

The explanation of whistler propagation is apparently 
successful in explaining the observed dependence of polariza- 
tion on kx angle, but it remains to see whether other 
observed characteristics are consistent with this interpretation. 
Clearly, the solar wind velocity, density, field magnitude 
and propagation direction should determine the observed 
frequencies if the waves are whistlers, so it is of interest 
to determine how sensitive the frequency is to variations in 
these parameters. 

Figure 12 is a figure similar to Figure 11, only it is 
plotted for the more limited frequency range of interest 
and for two sets of observed parameters. The solid curves 
correspond to density values of 3/cc and the dashed curves 
to densities of 12/cc. The lower two sets of curves represent 
the dispersion relations with the upper curve in each set 
corresponding to a field propagation angle of 20° and the 
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lower curve to an angle of 40°. The upper four sets of 

curves correspond to phase and group velocities calculated 

from the dispersion curves at the bottom with the two curves 

in each set corresponding to the two propagation directions. 

The vertical lines denote the regions where V , < 400 km < V 

ph g 

for each set of parameters. If the waves are to continually 
satisfy this condition (and there is no "a priori" reason 
they should do so) they must assume lower frequencies (and 
longer wavelength) for low densities and higher frequencies 
(and shorter wavelengths for high densities. The curves of 
Figure 12 are all drawn for a field strength of 5y but it is 
easily shown that for a given k the plasma frame frequency 
is proportional to the field strength. This means that such 
curves calculated for individual cases while always approach- 
ing zero at low frequencies will increase or decrease in 
proportion to the field strength. This fact plus the relative 
sensitivity of the curves to density and propagation direction 
make it quite plausible that the scatter of points in Figure 
8 is caused by variations in these parameters. 
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Another observed characteristic of the upstream wave is 
the fact that left handed waves were not found to exceed 

2 hz, whereas right handed waves were seen at least as high 
as 4.5 hz and perhaps extending to higher frequencies. This 
observation can be explained in terms of the whistler inter- 
pretation when it is recalled that the left hand frequency is 
produced by the V/\ term in equation 2'. Since for typical 
solar wind conditions we found that X was near 100 km, this 
term is typically of the order of 4hz. It in turn is 
decreased by adding the first term of 2' which is of the 
opposite sign, and typically has a magnitude between 2 and 

3 hz. Thus the sum is typically between one and two hertz 
and only for unusually high solar wind velocities or unusually 
short wavelength waves can the left handed waves appear at 
higher frequencies. In fact all of the left handed cases in 
figure 8 correspond to rather average solar wind densities 
and field strengths and therefore would not be expected to 
have higher frequencies. The highest velocities associated 
with left handed waves were 500 and 660 km/sec and these 

i 

correspond to 2 of the 3 highest frequencies (most 
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negative f's) observed near the shock in Figure 8. 

It is also of interest to look more closely at the 
one point in Figure 8 that is in apparent disagreement at 
f = 1.1, kx = 16°. This point corresponds to an unusually 
low density of 3, a field strength of 4.3y and a field- 
propagation angle of 49°. The simul.tanious solution of 
equations 1 and 2 reveals that the observed frequency of 
1.1 hz is the sum of an unusually large plasma frame 
frequency of 10.0 hz and a doppler shift term of -8.9 hz 
whose large magnitude is due to an unusually short wave- 
length of 50 km. The corresponding phase velocity is an 
unusually high 520 km/sec which is larger than the simultaneously 
observed solar wind velocity of 470 km, thus explaining the 
absence of a polarization reversal and why the right handed 
polarization was observed. 

It should be pointed out that the whistler propagation 
explanation for the observed waves is not dependent on the 
shock other than through the assumption that the waves are 
propagating away from the shock. It is also true that little 
has been said about why bands of frequencies are observed 
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other than to note that the plasma frame frequencies are 

near the lower hybrid frequency but not necessarily at it. 

One possible explanation of preferred frequencies is that 

they are generated with wavelengths which correspond to the 

scale of the shock thickness (e.g., Tidman and Krall, 1971). 

Full consideration of this question will be deferred for 

future work, but here we note that theoretical estimates 

2ttc 2ttc 

of the thickness of the shock or are of the order 

0) 0) . 
e l 

of 15 km and 700 km and thus span the range of wavelengths 
for the observed whistler waves. 

To further check the plausibility of this idea the 
simultaneous solution to equations 1 and 2 have been obtained 
for three interplanetary shocks where unusually narrow band 
upstream waves can be studied. Two of these shocks have 
been studied by Chao et al . (1973) and their orientation 
and velocity in the solar wind frame are known. Once the 
velocity of the shock relative to the spacecraft is known, 
a thickness can be calculated by noting the traversal time 
as deduced from the magnetic field data. This shock thickness 
can be compared with the independent calculation of the 
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wavelength. The relevant parameters are listed in table 1 
where At is the duration of the crossing and L the shock 
thickness . 


Table 1. Interplanetary Shocks 


Date 


At 

(sec) 

L (km) 

2nc <^ 

*e 

2r,c rkn,i 
U)^ f(hz) 

f'(hz) X (km) 

April 21, 

16:22 

3.0 

1380* 

12 

507 

1.0 

.5 

240 

May 17, 

6:25 

<.08 

<41 

8 

325 

4.2 

2.3 

89 

May 30, 

7:33 

.36 

162 

10 

442 

1.9 

.5 

170 


♦The unknown shock speed has been estimated as 
the solar wind speed plus 100 km/sec 


There is a general correspondence between X and the shock thick- 
ness with the smallest X corresponding to the smallest L and 
the largest X corresponding to the largest L. 

As another check on the whistler theory, we note that waves 
propagating away from interplanetary shocks will have a component 
of their k vector in the solar wind direction except perhaps 
in cases where the outward shock normal makes a large angle relative 
to the radial from the sun. Having k and V in the same 
general direction is the opposite situation from that of the 
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bow shock and requires that the sign of the last term in 
equation 2' be reversed. This implies that only right hand 
polarization would be expected to be observed ahead of inter- 
planetary shocks with conventional geometry. All three measure 
ments of the polarizations ahead of interplanetary shocks 
have in fact been right handed. 

Another question concerning upstream whistler waves is 
that of their damping. Figure 4 demonstrates that damping 
often takes place within a few minutes of the shock or, 
alternatively, in something like 100 cycles of the wave. 

There appears to be a tendency for larger field-normal angles 
to correspond to higher frequencies so one possible explanation 
of the slope in figure 4 is that the damping generally takes 
place over a given number of cycles but for larger Fn angles 
the cycles are shorter and the time interval during which the 
waves persist is smaller. This cannot be the complete explana- 
tion, however, because instances of high frequency waves 
occurring near spacecraft apogee in the sunward direction 
are found on occasion. In these instances, if we assume the 
waves are the same as those found near the bow shock, they 
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must have propagated many earth radii and undergone minimal 
damping over many hundreds of cycles . 

Cyclotron and Landau damping are candidates for damping 
upstream whistler waves since the propagation is not field 
aligned and the electric field will have a component both 
parallel and perpendicular to the background field. In these 
resonant interactions between waves and electrons, energy can 
be exchanged between the waves and the particles. Extraction 
of wave energy from the part icles results in damping but it 
should be noted that the energy of electrons moving away from 
the bow shock can also be converted into waves under the 
appropriate circumstances. In this sense it is not completely 
clear that whistler mode waves are propagating from the bow 
shock rather than being generated by electrons moving upstream 
from the shock. 

The question of non-linear Landau damping of electrostatic 
waves has been studied by O'Neil (1965) and the non-linear damp- 
ing of whistler waves has been studied by Palmadesso and Schmidt 
(1971) and Palmadesso (1972) . These results offer a promising 
explanation of the lack of damping. These authors have shown 
that after an initial partial damping of a large amplitude 
waves, the electron distribution becomes modified in a 
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manner such that it is unable to easily extract additional 
energy from the waves. For a monochromatic wave the damping 
ceases completely and for a more realistic wave spectrum the 
damping times are significantly lengthened over that expected 
on the basis of linear theory. If these theories prove applic- 
able it is perhaps more appropriate to ask the question of why 
the waves are usually damped quickly rather than persisting for 
long distances. 

Discussion . In reexamining the previously reported results 
in the light of the current observations it is clear that the 
waves discussed in the present paper are of the same type as 
those first noted by Heppner et al . (1967). These authors 
distinguished between one class of waves which occurred in 
packets and exhibited a rather sharp cutoff in frequency between 
.5 and 2.5 hz and a second class of more continuous waves at 
higher frequencies. The waves with frequencies less than 2.5 hz 
can be identified with the left handed waves of the present study 
and the higher frequency waves with the right handed waves. 
Gradations between the two categories which were not evident in 
the earlier study are now apparent so there no longer appears 
to be good reason to distinguish between these two 

categories, particularly since the whistler explanation seems applic- 
able to all waves. It is true, however, that left handed waves 
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often, but not always, occur in packets. This fact may be re- 
lated to the fact that such wave packets invariably occur in the 
presence of low frequency waves which are continuously 
changing the orientation of the "ambient" field. It also may 
be due to the fact that whistler waves are acting in the shock 
formation process in the manner described by Hasegawa (1972) . 

The packets away from the shock studied by Russell et al. 
(1971) may be a different phenomena than the waves studied 
in the present paper. The theory of Wu (1972) is more applic- 
able to these packets, whereas the theory of Hasegawa (1972) 
appears to be more applicable to the waves of the present 
study. The three cases of Holzer et al . (1972) exhibited 
right hand polarization and occured with the upstream field 
nearly perpendicular to the solar wind direction and hence 
they fit naturally into the results of the present paper. 

IV. Summary and Conclusions 

Waves of frequencies 0.5 - 4. hz persist upstream of 
the earth’s bow shock for intervals which are usually less 
than a few minutes in time but which occasionally are much 
longer. The waves propagate obliquely to the field at an 
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angle which is typically between 20° and 40 . The waves ,? 
may exhibit either left handed or right handed polarization 
with the particular polarization being controlled by the , 
propagation direction (alternatively the field direction) 
relative to the solar wind direction. The power spe.ctra of 
left handed waves decreases sharply beyond about 2 hz and 
predominant frequencies of left hand waves have never been 
seen beyond 2.5hz. Spectra of right handed waves often have 
peaks below 4 hz, but the spectrum decreases less rapidly 
than left handed waves. Similar waves appear to be present 
ahead of interplanetary shock waves. 

The cold plasma dispersion relation for right handed 
electron whistler waves is shown to be consistent with the 
observations in all respects. This equation predicts that 
waves of the observed frequencies could be detected upstream 
if they were emitted from the bow shock. It further predicts 
that such waves should be either left or right hand polarized 
depending on the propagation direction relative to the solar 
wind in just the manner observed. The left handed waves 

i 

are due to doppler shifting between the plasma frame of 
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reference and the observation frame. The 2.5 hz upper limit 
for left hand waves is due to the fact that wavelengths 
are invariably near 100 km for normal solar wind conditions, 
whereas unusually small wavelengths or large solar wind 
velocities would be needed to produce greater doppler shifts 
and higher left handed frequencies. It is tentatively suggested 
that the observed spectral peaks correspond to waves whose 
wavelength of approximately 100 km is related to the shock 
thickness . 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


- Illustrating the two classes of waves often seen 
upstream of the earth's bow shock. 

- Illustrating a "clean" shock crossing without 

low frequency waves (top) and a "pulsation" cross- 
ing in the presence of low frequency waves. 

- The number of "clean" and "pulsation" shock 
crossings as a function of the angle between 

the field and the shock normal. "Clean” crossings 
are confined largely to angles greater than 50° 
and pulsation crossings and upstream waves to 
angles less than 50°. 

- Persistence time of high frequency waves adjacent 

to the bow shock as a function of the angle between 
the field and the shock normal. 

; 

- Power spectra of a field component transverse to 
the average field direction for the examples of 
Figure 2. The April 12 example exhibits left 

.* ' t , 

handed polarization whereas the April 20 case exhibits 
right handed polarization. 



Figure 6 - Observed frequency of upstream waves as a function 


of the angle between the field and the approximate 
solar wind direction. Positive frequencies 
correspond to right hand polarizations and negative 
frequencies to left hand polarization. 

Figure 7 - Number of cases of upstream waves propagating at 
various angles to the upstream field. 

Figure 8 - Observed frequency of upstream waves as a 

function of the propagation direction relative 
to the approximate solar wind direction. Positive 
frequencies correspond to right hand polarization 
and negative frequencies to left hand polarization. 

Figure 9 - Clean shock crossing on April 3, 1971. The 

data appear to be typical but in fact are unusual 
since the upstream waves are left hand polarized 
rather than right hand polarized as in most clean 
shock crossings. 

Figure 10- Illustrating the presence of waves ahead of two 

interplanetary shocks. Calculated shock thickness 
and wavelengths are 162 km and 170 km on May 30 
and < 41 km and 89 km on May 17. 

yc 



Figure 11 - The cold plasma dispersion relation (solid curye'i 
and associated phase and group velocities (dashed 


curves) for typical solar wind conditions. The 
shaded region indicates where left and right 
handed whistler waves could be observed upstream 
from the earth’ss bow shock. 

Figure 12 - Cold plasma dispersion relations (bottom) and 
phase and group velocities (top) for two values 
of upstream densities. The shaded areas represent 
the regions bounded by field-propagation angles 
of 20° (top line of each pair) and 40°. Vertical 
lines help designate the regions where group 
values are greater than a typical solar wind 
velocity and phase velocities are less than the 
solar wind velocity. 
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